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Abstract. We consider atoms interacting each other through the topological 
structure of a complex network and investigate lattice vibrations of the system, the 
quanta of which we call netons for convenience. The density of neton levels, obtained 
numerically, reveals that unlike a local regular lattice, the system develops a gap of 
a finite width, manifesting extreme rigidity of the network structure at low energies. 
Two different network models, the small-world network and the scale-free network, 
are compared: The characteristic structure of the former is described by an additional 
peak in the level density whereas a power-law tail is observed in the latter, indicating 
excitability of netons at arbitrarily high energies. The gap width is also found to vanish 
in the small- world network when the connection range r = 1. 
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A variety of network systems such as computer networks, neuronal networks, 
biochemical networks and social networks, possesses complex topological structure, 
which can be described neither by regular networks nor by completely random 
networks |T|. Among various models for generating such complex networks, the Watts- 
Strogatz (WS) model j2] and the Barabasi- Albert (BA) model P provide the two most 
representative ones. Both model networks are characterized by very small network 
diameters proportional to the logarithm of the network size, which is dubbed "small- 
world behavior" and commonly observed in various real networks. Other important 
properties of many real network systems include high clustering and scale-free degree 
distributions jSj. While most existing studies of complex networks have focused on the 
geometrical and topological characterization, there now emerges an increasing number 
of studies paying attention to the dynamics defined on the networks. 

Lattice vibrations and associated phonon excitations on local regular networks, e.g., 
one-dimensional (ID) chains or two-dimensional square lattices, have been a textbook 
example of the introductory solid-state physics. In the case of such regular networks, 
atoms located on vertices interact only with nearest neighboring atoms, resulting in the 
standard linear phonon dispersion. On the other hand, when the interaction between 
atoms takes the form of the connection topology of a complex network, ingredients of 
long-range nature enter into the system. For example, the WS network has both local 
edges and long-range shortcuts j2] , which in this work are interpreted as the coexistence 
of the local and the long-range interactions. Existing studies on phonon excitations have 
mostly been performed on the local regular networks such as dimensional hypercubic 
lattices or on the fractal geometry jlj. Recently, the spectral properties of lattice 
vibrations have been studied in a small- world network similar to the WS network [Sj. 
However, effects of the local interaction range and the dependence of the spectral 
properties on the rewiring probability have not been investigated in detail. 

The main purpose of this work is to study in detail the vibration spectra of both 
the WS network and the BA network, in comparison with those of regular lattices. 
We thus consider atoms located on vertices of a complex network and interacting with 
each other connected via edges of the network, and investigate lattice vibrations and 
phonon excitations, with emphasis on the effects of the network topology. The results of 
such study can have implications in the mechanical property of a system of long chains 
bent in a complicated way. For example, in a bundle of long flexible polymer chains, 
some monomers which were separated by a long distance along a chain can make new 
couplings, building shortcuts. Throughout this work, we refer to phonons in a complex 
network as netons, to manifest qualitatively different characteristics due to the fact 
that the underlying geometric structure is not a usual periodic lattice. By means of 
numerical diagonalization of the dynamic matrix, neton excitation spectra are obtained 
and the corresponding density of neton levels are computed. It is revealed that the 
complex network in general has a finite gap in the level density, indicating the absence 
of low-energy excitations. The dependence of the gap width on the rewiring probability 
and the local interaction range is investigated in detail. The existence of the gap is 
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also reflected by the vanisliingly small specific heat at low temperatures. In the WS 
network three singularities are observed in the neton level density: Among them two 
are shown to be reminiscence of the van Hove singularity present in the local regular 
network, while the third one originates from the complex network structure. The BA 
network is also considered and the scale-freeness results in a power-law distribution of 
the level density, suggesting excitability of netons with arbitrarily high energies. 

We build the WS network following Ref. [2]: First, a ID regular network with 
only local connections of range r is constructed under the periodic boundary conditions. 
Next, each local edge (or link) is visited once, and with the rewiring probability P 
removed and reconnected to a randomly chosen vertex. After a whole sweep over the 
entire network, the average number of shortcuts in the system of size is given by 
NPr. In the WS network built as above, an atom is put on every vertex whereas an 
edge connecting two vertices represents the coupling between the two atoms located at 
the two vertices. For simplicity, we assume that all atoms are identical, each having 
mass M and moving only along the direction of the chain. The equation of motion for 
the Ith atom at the position xi then reads 

Mxi = Cj2i^rn-Xl), (1) 

meA-i 

where A; stands for the set of vertices connected to the Ith vertex (via either local edges 
or shortcuts) and the interaction strength C is assumed to be the same for any pairs of 
the interacting atoms. 

In the absence of shortcuts, Eq. describes a ID chain of atoms with only local 
connections up to the rth nearest neighbors. The analytic solution in this limit of P = 
is easily found to take the form: x/ oc e*('^'°~'^*), where a represents the interparticle 
spacing, u is the eigenfrequency and the periodic boundary conditions, xn+i = xi, 
restrict the wavenumber k to take the value k = 27in/N with integer n. For simplicity, 
we set a = 1, and obtain 



COS jk), (2) 



^ M 



which leads to the dispersion relation 



\ QM ^ ^ 

for small wave numbers (A; — > 0). Equation Q shows that u cc k regardless of the local 
interaction range r. Such linear dependence of the frequency u on the wavenumber k 
implies that the group velocity equals the phase velocity, both of which are frequency 
independent. Henceforth we measure the frequency u in units of v^C/M for convenience. 

In the presence of shortcuts (P ^ 0), the plane wave solution fails, owing to the 
lack of the translational symmetry. Substituting xi = xie^'^*' in Eq. (jT]), we write the 
equation of motion in the matrix form 

a;2x = DX, (4) 
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Figure 1. Frequency spectra of the WS network with the size N = 400 and the 
local interaction range (a) r = 1; (b) r — 2. The neton frequency uj is plotted as a 
function of the quasi-wavenumber q = 2Tm/N , where n (> 0) is the integer labelling 
the eigenfrequency from the smallest one. In the absence of shortcuts {P — 0), the 
linear relation ut r-^ q is shown for both r — 1 and r = 2. In the presence of shortcuts 
(P 7^ 0) , disclosed are very distinct frequency spectra between the two cases r = 1 and 
r = 2. 



where X is the A^-dimensional column vector with the components Xi {i = 1, ■ ■ ■ , N). 
The N X N dynamic matrix D has the elements 

{|Ai| forz = j, 
-1 forjGAi, (5) 
otherwise, 

where |Aj| denotes the number of vertices connected directly to the ith vertex via local 
edges or shortcuts. It should be noted that Dtj lacks the translational symmetry due to 
the existence of shortcuts. We have performed numerical diagonalization of the dynamic 
matrix D for various network sizes = 100, 200, 400, and 800, only to find insignificant 
difference between the two largest sizes. The eigenfrequencies found numerically for 
given network realization are sorted in the ascending order and an integer label n is 
attached from the lowest eigenfrequency {n = for the smallest u, n = 1 for the next 
smallest u, and so on). To obtain better statistics, we have also considered 1000-10000 
different network realizations, over which averages are taken. 

Figure H shows the neton frequency u versus the quasi-wavenumber q = 2'Kn/N at 
various values of the rewiring probability P for the local interaction range (a) r = 1 and 
(b) r = 2. Note that the quasi-wavenumber q becomes identical to the true wavenumber 
k only when P = 0. In general, the wavenumber k is not appropriate for the system in 
the presence of shortcuts (P 7^ 0): It does not make a good quantum number for P 7^ 0, 
where the translational symmetry is broken. For the range r = 1, Fig.^^a) suggests that 
the linear dispersion relation ~ g, which has been obtained analytically in Eq. for 
P = 0, still appears to hold at any value of P. In contrast, for the local interaction range 
r = 2, the dispersion uj{q) changes its functional form as P is raised from zero, reducing 
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to the linear relation a; ~ g only in the limit of P [see Fig. C^b)]. In particular, 
it is evident in Fig. ^b) that for P 7^ there appears a nonzero minimum value of u, 
below which the eigenfrequency does not exist. We have also examined the cases with 
the range larger than two (r > 2) and confirmed that these features remain qualitatively 
the same. With regard to this, it is of much interest to note that the WS network has 
a finite clustering coefficient ^ , only when r > 2 |2j , which makes it plausible that the 
difference in the neton frequency spectrum between the cases r = 1 and r = 2, shown 
in Fig. [TJ^a) and (b), respectively, is a manifestation of the clustering effects. We have 
also considered numerically the decomposition property of the network and found that 
while the WS network with r > 2 does not decompose into distinct parts, for r = 1 the 
probability of decomposition is finite. However, it should be noted that those networks 
with r = 1 which are not decomposed still display the vanishing gap. Consequently, the 
vanishing gap in Fig. ^a) for r = 1 has its origin not in the decomposition property but 
presumably in the clustering property of the network. Similar WS networks can also be 
built in a slightly different manner: Instead of rewiring local connections with the total 
number of edges fixed, one may add shortcuts to the network without removing local 
connections. In this case, we have confirmed numerically that there exists a gap even for 
r = 1, which keeps parallel to the observation that the XY model on the WS network 
with r = 1, constructed in such a way, has a finite-temperature transition, similar to 
that in the original WS network with r > 1 . 

We next consider the neton-level density g{uj), which is given by [S] 

g(ij) (X J dk 5{tu — uj{k)) (6) 

in one dimension. Using the mathematical identity 6{uj — uj{k))= \duj /dk\^^5{k — k^) 
with a;(fco) = we write Eq. (jH)) in the form 

(7) 

feo 

In the regular network (P = 0), Eq. (|7j), together with Eq. 0, gives the density g{uj) 
independent of cj, i.e., (7(0;) ~ tu^, in the low-frequency limit — > 0. 

In the presence of shortcuts, we use the box counting method to compute the 
neton-level density from the obtained eigenvalues of the dynamic matrix. The resulting 
density g{oj) is plotted in Fig. |2l which exhibits several characteristic features. First, as 
already found in Fig. d^b), the density g{uj) develops a gap ojc below which g{uj) = 0. 
Figure El shows that the gap uJc, well described by the linear relation ujc (x P for small 
P (see the inset of Fig. E}, keeps increasing as the rewiring probability P is raised and 
then appears to saturate for P > 0.7. Such saturation behavior has also been observed 
in the synchronization phenomena in the WS network [Oj. The emergence of the gap 
implies that the WS network is very rigid against lattice vibrations at small energies. 
In other words, one cannot excite netons at sufficiently low energies. We also compute 
the specific heat of the WS network for P = 0.1, and compare it with that of the local 
regular network (P = 0) in Fig. HI As expected, the presence of a finite gap in the 
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Figure 2. Ncton-level density g{uj) of the WS network with size N = 400 and the local 
interaction range r = 2 as a function of the frequency lu and the rewiring probability 
P. 



3 




0.2 0.4 0.6 0.8 

P 



Figure 3. Dependence of the gap width ojc on the rewiring probabiUty P. The inset 
demonstrates that Wc is well described by the relation Uc oc P for small P. 




Figure 4. Specific heat Cv of the WS network with A'' = 400 and r = 2 at P = and 
0.1. The existence of a gap in the neton- level density at P = 0.1 is reflected by the 
exponentially small specific heat at low temperatures. 
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Figure 5. (a) Neton-level density ^(aj) in the WS network with size N — 400 for r — 2 
and P = 0.01. (b) Heights of the three peaks in (a) versus P. Among the three, two 
peaks are the reminiscence of the local regular network and consequently their heights 
decrease with P. In contrast, the height of the third peak increases with P, suggesting 
that it has its origin in the complex network structure. 



neton-level density gives rise to exponentially small specific heat at low temperatures. 
This is in contrast to the local regular network, where thermal excitations of phonons 
are gapless and accordingly, the specific heat is described by the well-known Debye form 
with power-law dependence on the temperature, i.e., c„ oc T in one dimension. 

The second interesting feature observed in the neton-level density in Fig. |21 is the 
emergence of sharp peaks at small P. Figure Efa) shows in detail the density g{oo) 
at P = 0.01, displaying three peaks. It is well known that a local regular network 
(P = 0) (i.e., a regular lattice) has van Hove singularities whenever du/dk vanishes [see 
Eq. (|7j)]. For example, in the local regular network with the second-nearest neighbor 
interaction (r = 2), the positions of the van Hove singularities are easily obtained: 
The eigenfrequency uj = a/2(2 — cos k — cos 2k) given by Eq. combined with the 
condition du/dk = 0, leads to two singularities a.t u = 2 and u = 5/2; this allows us to 
conclude that the first two peaks in Fig. Efa) are remnants of the van Hove singularities 
present in the local regular network. On the other hand, the third peak at a; ~ 2.6 
cannot be explained in terms of the singularity of the local regular network. In order to 
probe the origins of the peaks, we measure how the peak heights change as the rewiring 
probability P is varied and plot the results in Fig. Efb). The first peak and the second 
one at u; = 2 and a.t u = 2.5, respectively, are observed to have heights which decrease 
with P. This again confirms that these two peaks originate from the local regular 
network, thus becoming weaker as the network loses its local structure. The height of 
the third peak, in contrast, is found to grow as P is increased, which clearly indicates 
that the third peak provides a feature generic to the WS network structure. 

Finally, we examine neton excitations on the BA network [Hj. In more detail, we 
first grow the BA network, following Ref . [H] : Starting from M vertices initially (at time 
t = 0), we attach a vertex with M edges at every time step. Once the BA network 
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Figure 6. Neton-level density g{u!) in the BA network with size N — 200 and AI — 3 
is compared with that in the WS network with N = 400, r = 2, and P — 0.1. The 
log-log plot of g{uj) versus oj discloses the power-law tail at high frequencies in the BA 
network. 

is built, we then put an atom on each vertex, assuming that the network is embedded 
in a ID chain. The frequency spectrum is obtained in the same way as that of the 
WS network, with the dynamic matrix still given by Eqs. Q and (0). Similarly to the 
WS network, revealed is the emergence of a gap. As the number M becomes larger, 
the gap grows and the peaks tend to smear out. Another interesting observation in 
the BA network is the slow decay of the density g{uj) at large frequencies, which is 
to be compared with the vanishingly small density at high frequencies found in the 
WS network. The log-log plot of the density g{uj) versus the frequency u, displayed in 
Fig. ini reveals that the BA network exhibits power-law decay which is well scaled as 
g{uj) ~ uj^'^'^. In contrast, g{uj) in the WS network shows very rapid decay, which is 
presumably of the exponential form. The observed power-law decay of g{uj) in the BA 
network implies that netons may be excited with arbitrarily high energies. It is of interest 
to compare this with the behavior of the adjacency matrix of the BA network 
Although the adjacency matrix A, with = and Aij = —Dij for i ^ j, and the 
dynamic matrix D have different diagonal components, both show similar scale-free 
distributions at large eigenvalues. 

In summary, we have investigated the phonon spectra in the monatomic chain where 
atoms interact with each other through the edges of complex networks. The density of 
phonon levels in the complex network, i.e., neton levels has been found to develop 
a gap, which indicates that netons may not be excited up to the critical frequency, 
and accordingly, a finite amount of energy is required to generate lattice vibrations. 
This robustness of a complex network has also been shown to yield exponentially small 
specific heat at low temperatures, in contrast to the power-law temperature dependence 
in the local regular network. The characteristic structure of the WS network has been 
described by the emergence of a peak in the level density in addition to those originating 
from the van Hove singularities. In the BA network, on the other hand, the level density 
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has been found to decay algebraically at high frequencies, which suggests excitability of 
arbitrarily high-energy netons. 
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